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h— -} • Abstract 

' We introduce the regularization Levi-Civita parameter for the 'curved Kepler', i.e., motion 

under the 'Kepler-Coulomb' potential in a configuration space with any constant curvature and 
metric of any signature type. Consistent use of this parameter allows to solve the problem of 
motion (orbit shape and time evolution along the orbit), thereby extending the use of the Levi- 
Q Civita parameter beyond the usual Kepler problem in a flat Euclidean configuration space. A 

'universal' description, where all relations are applicable to the motions in any space and with 
' any energy follow from our approach, with no need to discuss separately the cases where the 

, configuration space is flat or where energy vanishes. 

We also discuss the connection of this 'curved Kepler' problem with a geodesic flow. The 
i i ' well known results by Moser, Osipov and Belbruno are shown to hold essentially unchanged 

beyond the flat Euclidean configuration space. 'Curved' Kepler motions with a fixed value 
of the constant of motion —(2E — KiK2v7 2 ) on any curved configuration space with constant 
curvature k\ and metric of signature type K2 can be identified with the geodesic flow on a space 
with curvature a and metric of the same signature type. 

OO ■ 

^ ■ 1 Introduction 

l> ; 

For the Kepler motion in the Euclidean plane, the simplest way to integrate the Newton's equations 
involves an apparently artificial construct: by replacing the physical time t by a particular fictitious 
time s, called Levi-Civita regularization parameter, closed (and rather simple) explicit expressions 
for the cartesian x, y and polar coordinates r, <f> and, more importantly, for the time t can be given 
as functions of s (see e.g., Milnor [23]). The dependence in s of all these functions is smooth, so 
;_h ' this parameter in addition provide a way to regularize the description when the orbit approaches 

a collision orbit (where position and velocity are linearly dependent vectors) . The same parameter 
s appears also very directly in the relation between the Kepler problem and the geodesic flow on 
spaces of constant curvature [2H [23 [2S1 El HI H31 H] , and it turns out to be proportional also to 
the 'classical' parameter for studying Kepler evolution, called the eccentric anomaly (for elliptic 
orbits this parameter was introduced by Kepler himself [17]). Hence the Levi-Civita regularization 
parameter plays a central role in the Kepler problem. 

Our aim in this paper is to find the regularization parameter for the 'curved' Kepler problem, 
this is, the analogue of the Kepler problem when the configuration space is no longer the Euclidean 
space, but has a (non-zero) constant curvature. We perform this task using a Cayley-Klein type 
description, which allows to deal in a single run with the Kepler problem in a space <Swre 2 ] depending 
on two real parameters Ki and k%. The constant curvature of this space is K\ and its metric is either 
Riemannian, degenerate or Lorentzian, corresponding to the three alternatives K2 >,=,<0 for the 
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second parameter K2, referred to as the 'signature type'. As an added bonus of this approach, a 
complete and fully explicit solution for the 'Kepler problem' in a Lorentzian configuration space 
(DeSitter or Minkowski) is obtained. 

There are some papers dealing with particular aspects of either the classical or the quantum 
Kepler problem in configuration spaces of constant curvature but none of them (to our knowledge) 
deals with the regularization parameter for the problem, nor its possible relation with a geodesic 
flow. And in spite of recent interest in studying motion in Lorentzian manifolds, we do not know 
cither any paper dealing with the 'Kepler' problem on a Lorentzian configuration space, curved or 
not. Thus, when the configuration space has nonzero curvature, or when the metric is of Lorentzian 
signature, the results obtained are new. 

We solve completely this problem, and recover as a particular instance all the well known results 
for the 'flat' Kepler problem in Euclidean space. In particular, the close relationship between the 
set of Euclidean Kepler motions with total energy E and the geodesic flow in a space of constant 
curvature — 2E [211 [25j [26J [H HI [23] , appears as the particular 'flat case' of a generic relation holding 
for the 'curved' Kepler motions in a configuration space S K r K2 i °^ an y cons t an t curvature and either 
metric signature type. If a denotes the combination a = —(2E — K1K2S 2 ) of energy and angular 
momentum, the result we obtain is: the set of 'curved' Kepler motions with a given constant value 
of er can be identified to a geodesic flow in a space of constant curvature a whose metric has the 
same signature type as the configuration space. This result extends directly the Euclidean Kepler 
one, yet we have not found any reference to this, nor to any similar result, in the literature. 

All the expressions we give are completely explicit; for the Kepler motion, we disclose the depen- 
dance of all the relevant coordinates on the 'curved' Levi-Civita parameter s. This automatically 
produces a large number of results and relations for the 'curved' Kepler problem, which extend prop- 
erties well known for the Euclidean Kepler problem: the 'cycle' character of the Kepler hodographs, 
the Kepler equation and the dependance t(s), etc. The connection of Kepler motions with fixed 
value of the quantity a with a geodesic flow in an auxiliar space whose curvature is precisely a is 
quite direct. 

This is done using as a language the parametric CK type approach, which allows to do compu- 
tations for all the configuration spaces at a single run, considering k\,K2 as free parameters. But 
there is more: even if we stick to studying the Kepler problem in, say, Euclidean space E 2 , where 
the two CK free parameters have fixed values n\ =0, K2 = 1, the quantity a, which is in this case the 
energy up to a factor, which plays a role as a third Cayley-Klein parameter, working consistently 
in these terms allows to give an unified description, for all types of orbits, in a single run, and with 
a dependence on a which is smooth when a — > 0. The same happens when the configuration space 
is the general £ rei r re2 ]- Thus the CK type approach provides some new perspective to an 'universal' 
formulation for the Keplerian orbits, encompassing not only all energies but also all possible values 
of the curvature and signature type of the configuration space. 

The plan of the paper is the following: Section 2 is devoted to the regularization of the 'curved' 
Kepler problem. First we give the basics on dynamics on a space S 2 ^ K2 ^ and introduce the 'curved' 
Kepler potential and the 'curved' Levi-Civita parameter s. Then we derive expressions for coordi- 
nates and time in the curved case as functions of s. Section 3 is devoted to the connection between 
Kepler motion in Sw^i and a geodesic flow. This is done directly, in terms of a stereographic 
projection, which identifies 'curved' Kepler motions with a fixed value of the constant of motion 
a = —(2E — K1K2J 2 ) to the geodesic flow on a space with curvature a and metric of the same 
signature type as the configuration space. Again this reduces, when the configuration space is Eu- 
clidean, to the well known Moser-Osipov-Belbruno result: Kepler motions in E 2 with energy E can 
be identified with the geodesic flow on a Riemannian space with constant curvature —2E. The three 
'curved' Kepler laws are stated and discussed in Subsection 3.4. Finally, Section 4 discusses the 
specialization to the Euclidean configuration space, translating the results to the standard language 
and relating the parameter s to the eccentric anomaly. 
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2 The regularization of the curved Kepler problem 



2.1 Dynamics in a configuration space S^r K2 , 

We denote by S K r K2 ^ a 2d space with constant curvature k% and metric of a signature type (1, k%). 
By rescaling lenghts and angles, the two parameters k\,K2 could be brought independently to a 
standard value 1,0,-1 and the nine combinations correspond to the so called 'standard' Cayley- 
Klein CK spaces. Instead of reducing m, K2 some standard values to start with, we will leave both 
Ki,K2 as free parameters. There are two main reasons to do so. First, a unique computation, 
only slightly more complicated than the one required for each individual instance, covers all nine 
cases. Second, dealing with general values for Ki,/t2, one can get a valuable perspective on how 
different properties and relations depend on the curvature and/or signature type of the space, and 
how things change when curvature vanishes or changes sign, or when the metric degenerates or 
changes from positive definite to indefinite. This makes a subsequent analysis of contractions and 
limiting cases completely redundant. For more details on this formalism, see [15| 15) IS]. 

The basic tool in this CK approach is the use of a set of 'labeled' trigonometric functions. 
The K-labeled 'Cosine' C K (x) and 'Sine' S K (x) functions are defined here as the solutions of the 
differential equation: 

dx 

determined respectively by the initial conditions: 

dS K (x) 



F(x) = -kF(x), (1) 



c K (o) = i, dC ^ x) 



dx 



= 0; S K (0) = 0, 



x =a dx 



= 1. (2) 

x=0 




In the field of abstract differential equations, a similar approach leads to the 'cosine' and 'sine' 
families [TT]. Here, with the simplest equation ([T]) this 'cosine' and 'sine' are ordinary functions, 
which admit the following expressions, with an analytic dependence in both variables x and k: 

-j= sin \/Tl x k > 

C, I.N := J i . S, (.,•] : { X _ K = . (3) 

—k= sinh J — k x k < 



These two functions satisfy the basic identity (as well as many others, see e.g. [15]): 

C 2 k (x)+kSI(x) = 1 (4) 

The 'Tangent' T K (x) is defined as the quotient T K (x) — S K (x)/ C K (x). Another function ap- 
pearing naturally is the k version of the 'versed sine', defined as V K (x) = (1 — C k (x))/k; it is 
interesting to realize that when n — > 0, both numerator and denominator tend to in such a way 
that V K (x) is well defined even when n — > and Yq(x) = x 2 /2. The inverse functions are denoted 
accordingly; here only ArcT K (ir) will appear. 

These functions allow us to write expressions in any S 2 ^ K2 ^ in a unified way. For the two 
standard choices k = ±1, these functions are precisely the circular or hyperbolic trigonometric 
functions. The singular case k = corresponds to the so-called parabolic trigonometric functions, 
and for general n, the labeled 'cosine' C K (x), 'sine' S K (x) and 'versed sine' V K (x) functions can 
be considered as an one-parameter set of deformations of the corresponding 'parabolic' functions, 
equal respectively to l,x and x 2 /2, to which these reduce for k = 0. Within this viewpoint, the 
tangent T K (x) is to be considered as another different deformation of the function x. 

The space S K1 i K2 \ involves two independent labels K\,K2, and in its CK description K-labeled 
functions with the labels «i, K2, K1K2 appear. When the two basic labels are positive, the space 
>S Kl r K2 i is a two-dimensional sphere; the standard S 2 corresponds to the choice Ki = l, K2 = l; other 
standard choices are k\ = 0,K2 = 1 (= E 2 , the Euclidean plane) or K\ = — 1,^2 — 1 (= H 2 , the 
hyperbolic or Lobachcwski plane). The remaining standard spaces are the three 'non-relativistic' 
(with absolute time) space-times, appearing for K2 — 0: antiNewton-Hooke space ANH 1+1 , Galilean 
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Table 1: The nine standard two-dimensional CK spaces S^_^ K2 y 





Measure of distance & 


Sign of ki 


Measure of angle 


Elliptic 


Parabolic 


Hyperbolic 


& Sign of K2 


K\ = 1 


Kl =0 


Kl = —1 




Elliptic 


Euclidean 


Hyperbolic 


Elliptic K2 = 1 


S 2 


E 2 


H 2 




Co-Euclidean 


Galilean 


Co-Minkowskian 




Oscillating NH 




Expanding NH 


Parabolic K2 = 


ANH 1+1 


qX+1 


nhw 




Co-Hyperbolic 


Minkowskian 


Doubly Hyperbolic 




Anti-de Sitter 




De Sitter 


Hyperbolic K2 = — 1 


AdS 1+1 




dS 1+1 



space G 1+1 , Newton- Hooke space NH , and the three relativistic space-times in 1 + 1 dimensions, 
appearing for K2 < 0: the AntiDeSitter sphere AdS 1+1 , Minkowskian space M 1+1 and deSitter 
sphere dS 1+1 ). This information is displayed in the Table; see for more comments. 

At a first contact with this formalism, the reader might well pretend that «i = 1,k% = 1, and 
all the labeled 'cosine' C K (x) and 'sine' functions S K (x) with either label (k±,K2 or k±K2, all equal 
to 1) are ordinary, circular cosines cos(a;) and sines sin(x). With this understanding, every relation 
will apply to the standard sphere S 2 , where it is possible to visualize most properties in an easy 
way. At the final intended level of reading, «i, K2 should be considered of course as free parameters, 
and the CK formalism keeps track automatically of all sign differences, vanishing of some terms, 
replacement of (some) circular trigonometric functions by their parabolic or hyperbolic analogues, 
etc., which distinguish relations in the nine CK spaces. 

As a first example of this language, let us write the expressions for the metric in the ^i[k 2 ] > m 
the intrinsic polar (r,<f>), parallel '1' (x,v) and parallel '2' (u,y) coordinates: 

dl 2 = dr 2 + k 2 S 2 Kl (r) ddp 2 = dx 2 + n 2 C 2 Kl (x) dv 2 = C 2 KiK2 (y)du 2 + k 2 dy 2 . (5) 

Comparison with the very well known form for the metrics in some standard spaces as S 2 , E 2 , H 2 or 
M 1+1 will give some feeling on the unification capacity of the CK formalism. In the general curved 
CK space the quantities u and x (or v and y) are generally different (this is clear on a sphere); this 
follows from the relations holding for the general CK space S k j K2 \ : 

T Kl (u) = T Kl (r)C K2 (4>), T KlK2 (v) = T M (r) S rt2 (0) (6) 

S Kl (x) = S Kl (r) C K2 {cf)), S KlK2 (y) = S K1 (r) S K2 ((j>) (7) 

which can be considered as formulas for orthogonal triangles |15j in the trigonometry of S 2 ^ K2 y 
When Ki = these relations imply a characteristic degeneracy of the flat spaces: the equalities 
u = x and v = y. From the point of view of deformations, © and ([7]) can be looked at as two 
possible 'curved' generalizations of the Euclidean relations x — r cos (f>, y — r sin <f>. Notice that 
(u, y) and (x, v) are orthogonal coordinates in any S^J^] > ( u ' u ) or ( x > V) are orthogonal systems 
only when the curvature k± vanishes; in this case all these four systems collapse to a single Cartesian 
system. 

Consider now the motion of a particle in the configuration space S 2 _^ K2 ^ under a natural mechan- 
ical type Lagrangian, with a kinetic term given by the metric and possibly a potential depending 
on the coordinates: 

C=^g^(q\q 2 )q^-V(q 1 ,q 2 ). (8) 

Constants of motion for this Lagrangian which are linear in the velocities occur only if the potential 
is invariant under some one-parameter subgroup of isometries, and are equal to the corresponding 
Noether momenta. With a scale factor suitable to simultaneously deal with all the spaces in the CK 
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family [5] , the natural base of three Noether momenta V\ , V2 , J is given in the three coordinates 
as: 

Vi\ ( C K2 (^)r - k 2 C Kl (r) S Kl (r) S K2 
V2 = S K2 (0)f + C Kl (r) S Kl (r) C K2 
^/ V S 2 Kl (r)0 

^1 \ / C KlK2 (v)x + kik 2 S KlK2 (v) S K1 (x) C Kl (x)v \ / Gl lK2 (y)ii 
V 2 = C* Kl 2 (a> = K X S Kl (u) S KlK2 (y) C KlK2 (y)u + C Kl ( U )y 

J/ \ -S KlK2 (v)± + C KlK2 (v)S Kl (x)C Kl (a;)?) / \ - C Kl (u) S KlK2 (y) C KlK2 {y)ii + S Kl (u)y 

In the standard Euclidean E 2 , where K\ = 0, K2 — 1 and both parallel type coordinates reduce to 
Cartesian ones, the CK momenta are in these coordinates: 

Vi\ w =± = u, V 2 \ ¥;2 =v = y, J\ E 2 = xv - vx = -yii + uy = xy - yx . (9) 
2.2 The curved Kepler problem 

The 'curved' Kepler potential in any S R < K i configuration space is: 

Vk = -=^- (10) 

In the three Riemannian spaces of constant curvature (CK spaces with K2 > 0), this potential 
follows by enforcing Gauss law in the corresponding CK three dimensional space, and in particular in 
the hyperbolic plane it was considered early by Lobachewski himself; see some historical comments 
in [10) . In quantum mechanics on the sphere case this potential was discussed in the paper by 
Schrodinger (27J , shortly to be followed by Infeld and Schild [TB] who studied the hyperbolic space 
case. 

Motion in Kepler potential in a curved configuration space has two trivial constants of motion, 
energy and the angular momentum: 

E = \{Vl + K2Vj + K l K2J 2 )--^— J = S 2 Kl {r)4>, (11) 

and there are two additional constants, to be discussed later, which are associated to the separability 
of the Kepler potential in two equiparabolic coordinate systems (see [?])■ 

In the , the orbit is known to be a conic (for the intrinsic geometry of S 2 _^ K2 ^) 7 with a focus at 
the potential origin [HUH]. In polar coordinates, with the periastron placed on the halfdine (j> — 
(hereafter taken as the standard position of any Kepler orbit), its equation in the CK space S K j K3 -\ 
is: 

T Kl (r) = V^T KlK2 (p) * (12) 
1 + e C K2 (</>) 

The two geometric parameters determining the conic, i.e., its CK eccentricity e and semilatus rectum 
p are related to the physical conserved quantities energy E and angular momentum J as: 

/ — rp ( , K2J 2 2 1 , (2E - Kl n 2 J 2 )n 2 J 2 
V«2 T KlK2 (p) = — — , e = H p (13) 

Orbits in general position involve a further constant </>o : the angular position of the periastron; its 
equation follows from (|12p by the trivial replacement 4> —></> — ipo . The constant </>o is not essential, 
as it comes from the central nature of the potential, and will be disregarded from now on. 
We note also a simple relation between r and u along any orbit: 

T Kl (r) = ^7— - e T K1 (u) = ^/k^T KiK2 (p) - e T Kl (u) (14) 

which follows from the relations © among coordinates and hold for any CK space; these lead to 
the ^-dependence of u, v for the orbit in standard position: 

T Kl (u) = V^T KlK2 (p) ^tttx , T KlK2 (v) = V^T KlK2 (p) ^ttttn ■ (15) 
l-eC K2 (0) l-eC« 2 (0) 
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2.3 Levi-Civita regularization of Kepler motion in a curved configuration 
space 

As the potential is central, the angular momentum constant J leads to the 'curved' version of the 
law of areas: 

In the euclidean case, time evolution x(t), y(t) for cartesian coordinates satisfies a non linear 
system of differential equations, which do not allow closed form integration in terms of elemen- 
tary functions. For elliptic orbits, the problem of motion was solved by Kepler [17j . essentially 
introducing a new parameter along the orbit, the auxiliar angle known as the eccentric anomaly, 
instead of taking the polar angle 0, which was known in the classical parlance as the true anomaly. 
The eccentric anomaly (which can be also suitably defined for parabolic and hyperbolic orbits) 
surprisingly simplifies the situation so as to allow closed elementary expressions. 

The simplification afforded in the Kepler problem by the use of the eccentric anomaly can be 
considered as a consequence of a more fundamental structure: when the time t is replaced by 
another parameter, the (fictitious) Levi-Civita time s (also known as the Levi-Civita regularization 
parameter) related to t by the condition s = l/r, the equations for the s-evolution of x,y become 
linear [20]. If this happens in the Euclidean Kepler problem, a natural question for the 'curved' 
Kepler problem is: does there exists a 'curved' regularization parameter s and some (functions) of 
coordinates on S k j K2 i (providing some curved extension of the cartesian coordinates x,y), whose 
s-evolution for the curved Kepler problem is given by a linear system? And to which extent and 
precisely how other properties already known for the Euclidean Kepler problem still hold for the 
'curved one'? This is the problem we address in this paper. 

Any ansatz for the 'curved' functions of coordinates and the 1 curved' regularization parameter 
should reduce as k\ — > to the cartesian coordinates and to the 'Euclidean' Levi-Civita parameter 
determined by s = l/r. 

Our choice for the curved version of the regularization Levi-Civita parameter is: 

. _ ds _ i r i 

S " dt ~ Sjrj Cjr) ' S ~ J S^Cjr) {U} 

At first sight, this seems to be an arbitrary choice over other apparently more natural possibilities, 
as s = 1/ T Kl (r) or s = 1/ S Kl (r); the results will justify however the particular choice made here. 
By combining with the law of areas we get as well: 



ds_ _ T Bl (r) 
deb ~ J 



s = iy"T K1 (r)#. (18) 



Direct computation leads then to the following result: Along any curved Kepler motion, the s 
evolution of the two coordinate functions T Kl (u), T KlK2 (r) © is given by the linear system 

-^-T Kl (u) = -^-T RlK2 (v) , ^-T KlK2 (u) = ej + -|- — —T Kl (u), (19) 
ds J as j K2 

hence satisfies the (uncoupled) linear second-order system 



d 2 T ( , , fc 2 l-e\ , . d* P 1-e 2 

— 2 T Kl (u) = -ke-——- T K1 (u) , — 2 T KlK2 (v) = -J5 — 



The coefficients of T Kl (u) , T KlK2 (v) in the r.h.s. are equal in the two equations and the sign 
of this quantity will determine the nature of the s-dependence for the solutions. Expressing this 
coefficient in terms of the energy and angular momentum through (|13j) . we obtain: 



J 2 K 2 



= -(2E - k iK2 J 2 ) . (21) 
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Here we see that in the curved case, the quantity governing the character of the evolution is not 
—IE as it was in the Euclidean case [53] but rather the combination — (2E — H1K2J 2 ), which has, 
further to the energy, a contribution from the angular momentum. If we define a as 

a := -(2E - Kl n 2 J 2 ) , (22) 

and arrange the constant terms, the second order system is: 

( T Kl (u) + ^)= -*(t ki (u) + ^) , ^ T KlK2 (v) = -a T KlK2 (v) , (23) 



whose general solution, for either ^T Kl (u) + or T KlK2 (v) as functions of s is, after (fTJ), a linear 

combination AC a (s) + BS a (s) of the two basic solutions C a (s), S CT (s). This is the more general 
solution, even when a = 0, for then it reduces to A + Bs. The symbol er has been chosen to 
underline its natural link with s, as the 'CK label' of the Levi-Civita parameter s (just as the 
lenghts u,r have 'label' K\ and the angles <f> have 'label' K2, in the sense any appearance of these 
quantities is through the corresponding 'labelled' functions). 

Now for the evolution along the orbit in the standard position (with </>o =0), at the periastron 
v must vanish and u must be extremal; if the origin of s is chosen also at the periastron, enforcing 
the previous conditions when s = 0, the expressions for T Kl (w), T KlK2 (w) are narrowed down to: 

T Kx (u) = A C a (s) - — , T K1K » = B S a (s) (24) 
a 

and finally A, B can be related to the physical constants along the Keper motion by noticing that at 
the periastron, where s = one must have u — r per , and hence by equating T Kl (u)| s=0 = A — ke/a 
to T Kl (r per ) = y/K2 T Klft2 (p)/(l + e) taken from (Ti"2"l) and using (|22|) we get after some algebra 
A = k/a. Hence: 

T K1 (u) = k Ca(s) ~ e . (25) 

(T 

To determine the constant B, let s s i be the value of the parameter s at the semilatus rectum of 
the orbit. At this point u — and v = p. We have then two equations: 

k C " (5 *' ) ~ e =0, B S CT ( Ssi ) = T K1K2 ( P ) (26) 
a 

and elimination of s s i through the basic identity C^(s) + crS^(s) = 1 leads, after some work to 
B = ±J. The choice of sign corresponds to the sense of motion along the orbit and in what follows 
we will choose B = J . Thus: 

T KlK2 (v) = JSa(s) (27) 

Up to now we have found the s dependence of u and v (through T Kl (w), T KlK2 (v)). The depen- 
dence of the radial coordinate r on s follows directly from 



T Kl (r) = - ke a(5) ~ 6 = k l - eG ^ (28) 

k a a 

where we used the relations (|13[) which are better rewritten in terms of a as 

V^T KlK2 (p) = fci-^. (l- e 3) = ^!. (29) 



a 



In both equations, when a — > or when K2 — > 0, then e — > 1 while the quotient (1 — e 2 )/(<jK2) = 
J 2 /k 2 remains well defined and depend only on the angular momentum Notice that in spaces with 
a degenerate metric {n2 — 0), all Kepler orbits have e = 1, in this case the orbit equation (fT2"| 
should be rewitten in terms of the versed sine of the angle <j). 

The use the Levi-Civita parameter enables the finding of closed exact and smooth expressions 
for the s-dependence of u, v and r in the general 'curved' case almost as easily in the Euclidean 
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case. For the relation among s and the angular coordinate <f> we may use (|18|) and replace T Kl (r) 
taken either from (fT2"|) or from (|28|1 to find two alternative forms: 

ds_ _ k_ 1 - e G a (s) _ J_ k 2 

dc/> ~ J a ~ k 1 + e C K2 (<?!>) 1 ' 

Each of these forms allows closed integration by using the formula: 

/ — -j— ds = - 2 ArcT_ CT ( f +1 T(j H 1 (31) 

7 l-eC CT (s) \v^t 

Integrating (|30[) (notice that the integration can be done in two ways as the equation allows two 
separated forms) and using (J2SJ) gives: 

^A=^T_J4±^Tji)) (32) 



-a 2 Vv^^T 
leading after some manipulations to the s-dependence of the angular coordinate 



^T K2 (|)=(l + e) ^T ff (|) (33) 

which remains meaningful even when a = 0, as clear from (|29p . This formula does not involve 
explicitly the curvature K\ of the configuration space. For the standard Euclidean case, where 
K2 = 1, the l.h.s. will be tan(</>/2), so this formula ressembles closely the well known Euclidean 
relation between the angle (f> and the eccentric anomaly; indeed there is a direct relation, in the 
Euclidean case, between s and the eccentric anomaly, as we will see in the last section. 

The equation (l3^|) displays a curious symmetry among C K2 (</>) and C CT (s), which can be enhanced 
by writing it in the form: 

l-eOOQ 1 + eCM = J^ = 1^ (34) 

G K2 k 2 (JK2 

Closed formulas can also be obtained for the three basic (cosine, sine and tangent) trigonometric 
functions of the angular coordinate; the most direct way is to recall the relations © which hold 
for the coordinates in any CK space obtaining: 



C K2 (</>) 



C CT (s) 



l-eC» 



o ,.x J crSgjs) / 1 - e 2 a Sgjs) 

°* 2 (W = t; Fryr = W i Fryr 35) 

k 1- e L a (s) V k 2 o 1 - e C^s) 



T K2 (0) 



J aScis) /l-e 2 aS CT (s) 



Co(s) - e U k 2 ct C CT (s) 



(of course, starting from these expressions, the relation (|33p can be also derived by direct compu- 
tation instead of by integration). 

We already have obtained the s-dependence of the radial coordinate. Let us ask now for the time 
t as a function of s. By starting from the definition (fTT|) , using the identity 1/ G 2 K (r) = l-\-K\ T 2 . (r) 
and the s-dependence of T Kl (r) we obtain: 



2 

ds 1 _ 1 _ l + Kl T 2 Kl (r) _ 1 + ^i^fc 



(36) 



dt S K1 (r)C Kl (r) T Kl (r)C 2 K1 (r) T Kl (r) k l^M 
so that, with the conventional choice for the time origin t = when the particle is at the periastron: 

ifl-eC^)) 



l + Kl ^(l-eC CT (s)) 



■ ds (37) 
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The integration can be done by decomposing in simple fractions, and leads to: 

1 f ( a — e ft + e 1 . . 

t[s) = J { T^aw + TTpcJs) } ds (38) 

where a, ft are dimensionless quantities: 



ke.\J— Ki ke^J—K 



ft = (39) 



a + kyJ—Ki a — k\J—K\ 

which by using (|3ip can be expressed in completely closed form: 

t(s) = _L= { ^== ArcT_ CT ( -J±L T a ( S -))+ 4=== ArcT_ ff ( J=== T a ( S -)) ) 

(40) 

By noticing a — e = — k J_ Ki a, f3 + e = k J 7 _ Ki ft, the previous expression can be also given as (notice 
the dissapearance of e in the previous a — e, ft + e terms, the relative minus sign and the change in 
the global prefactor in the new expression): 



ekm [v^^T ' Vv^^T \2JJ Vv 7 /^ - 1 

(41) 

Summing up: closed expressions for the coordinates u, v, r, <f> and the time t along a Kepler orbit 
in a curved configuration space have been obtained. Together with ([25l l27l |28|) . this expression 
completely solves the configuration space with any constant curvature and any signature, and for 
any value of the constant a. 



3 The Kepler motion in a curved configuration space as a 
geodesic flow 

Now we want to explore whether or not the well known relation among the Kepler motions with 
energy E in Euclidean space and a geodesic flow in a space of constant curvature —2E (and definite 
positive metric) extends for the 'curved' Kepler motion in any S%i K2 y A constructive way to 
establish this connection is to look to the Kepler evolution not in the configuration space S K < K2 \ 
but in an auxiliary space, the space of Cayley-Klein momenta Vi,V2- 

The 'curved' Kepler problem has, further to energy and angular momentum, two additional 
constants of motion, which can be considered as the components of a single vector, the eccentricity 
or Hamilton vector, which is related to the Laplace-Runge-Lenz vector (see [5] for comments on the 
relation between the eccentricity vector and the Laplace-Runge-Lenz vector in the 'curved' case and 
[121 H31 [21] for additional information in the Euclidean case). In any configuration space S^j Ka -i, 
[9] these constants are: 

£oi = JV\ + k S K2 (0) , £ 2 = JV 2 + k V K2 (0) (42) 

The existence of these Kepler first integrals can be seen as a consequence of the separability of 
the curved Kepler potential in two systems of 'parabolic' coordinates in the curved configuration 
space (with a focus at the potential origin) (7j. The particular form of the constants (|42|) has a 
geometric consequence: the 'momentum hodograph' is a 'cycle' in the plane of momenta Vi,V2, 
a plane whose metric has signature type Ki (this is suggested in Taking the constancy of 

(|42|) as the departure point would be a modern, natural choice. For the Euclidean Kepler problem, 
the circle character of hodographs was first stated explicitly by Hamilton, who actually derived 
things the other way round: from a suitable rewriting of Newton's equations, the circle character 
of the (velocity) hodograph follows, and this leads to the two new constants of motion specific to 
the Kepler problem, the Hamilton eccentricity vector. So let us mimic, for the general CK space, 
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the path trodden by Hamilton in the Euclidean Kepler problem [14] and let us try to derive the 
constancy of (|42jl . keeping in mind an important fact: Hamilton considered the velocity hodograph 
vector, and we are considering the 'momentum hodograph' T > \,T > 2\ this will be of some relevance 
later. 

When a particle moves in S K v K ^ following Kepler evolution, angular momentum J is constant, 
but V\,Vi do depend on time. Their time evolution is the extension to the space S Kl r K2 i of the 
original form of Newton's equations: the rate of change of momentum is equal to the force. This 
'force' is a vector under rotations around the origin (hence, a vector in a flat plane with signature 
type K2), and its components are (J 7 !,^) = — k(C K2 (<j)), S K2 (<f))/ S K1 (r), whose radial dependence 
1/ S„ (r) follows directly from the Gauss law in a 3d space of constant curvature m (where the area 
of the sphere of radius r grows as S 2 Kl (r); notice also d(k/ T Kl (r))/dr = -k/S 2 Kl {r))). Hence we may 
write: 

^-Tr-sS^ -**'- A = ^ = ~im s ' M - <43) 

Elimination of t using of the law of areas leads to: 

a system whose integration is trivial, leading precisely to (|4"!?j) . Had we started from the constancy 
of the two components of the eccentricity vector, the equations (f4~5]l and l|4"4*|) would follow simply 
by differentiating (j4"2|) with respect to time and using the law of areas (fT5]l . The orbit follows 
directly by enforcing the relation between the three CK momenta which generalizes the euclidean 
J = xV 2 -yVx UM- 

This means that the 'momentum hodograph' curve is always a 'cycle' in the Vi,V2 space rela- 
tively to the (flat) metric with signature type k 2 , based in the quadratic form dp 2 := dV\ + ^V 2 - 
(We recall cycles are defined as the curves with constant geodesic curvature; in the Euclidean plane 
these are circles, with straight lines as limiting cases). Of course in the Euclidean Kepler problem, 
the two momenta V\ , V2 are simply ^ equal to the two cartesian components of the ordinary 
velocity vector. Thus the Hamilton result for the Euclidean Kepler problem, formulated in terms 
of velocity hodographs, could equivalently be restated in terms of momentum hodographs without 
any change whatsoever. But this equality does not extend for curved spaces. No simple extension 
of the Hamilton result should be therefore expected for the 'velocity' vector in the 'curved' Kepler 
case because the velocity is a tangent vector to each point of the trajectory, which would require 
some kind of transport to a common origin to make sense of a 'curved velocity hodograph'; this 
would mean some further arbitrary choice as to the path to follow for performing the transport. 
On the other hand, the 'momentum hodograph' makes sense in any curved configuration space, 
without any further assumption, because the Cayley-Klein momentum Vx,V 2 i already lives in a 
Lie algebra and can be considered as a vector attached to the origin of a linear space without any 
need to perform any kind of 'transport'. 

Hence we may state as the first main result 

Theorem ['curved' version of Hamilton's] As t varies, the CK momentum vector V\,V2 of a 
particle undergoing Kepler evolution in a space S K1 r K2 :, moves along a 'cycle' in the 'momentum 
plane' whose metric dp 2 — dV\ + K 2 dV 2 IS an< ^ 0I ~ signature type «2- 

Now we ask about the relation among the Levi-Civita parameter s with the metric dp 2 in this 
result. The (square of the) 'norm' of V = dV /dt = (dVx/dt, dV 2 /dt) seen as a vector in a 2d plane 
whose (fiat) metric has signature K2 is: 

V -V = {V\f + = (45) 

Consider now the s-dependence of momenta, related to the time dependence as: 

dp dV dt dP 

-ds- = HdS = s - (r) c - (r) ~dt ' (46) 
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so that for the 'norm' of the vector dP/ds we have: 



* ^ - SKi(r)CKi(r) ^>)-^- (47) 

Conservation of energy pip allows to draw k/ T Kl (r) in terms of momenta, 

— k — =E- \h x k 2 J 2 - \ {{Vrf + k 2 (V 2 ) 2 } = -| {(Pi) 2 + n 2 (V 2 ? + a} (48) 

where again the quantity a appears. By combining these relations, we get: 

dg2 = 4{(^ 1 ) 2 + . 2 (^ 2 ) 2 } = 4 

((Pi) 2 + n 2 {V 2 y + a) 2 ({Vif + K 2 {V 2 f + a) 2 

Hence, if (|49|) is considered as the definition of a metric in the momentum plane, the increase of the 
Levi-Civita parameter along any Kepler orbit is equal to the 'lenght' provided by this 'Levi-Civita 
metric' (|4"5|) along the curve described by the evolution of the CK momenta V\,V 2 . The metric ds 2 
on the space of all CK momenta admissible for a Kepler motion with fixed a, does depend on a. It 
is different from the flat momentum plane metric dp 2 implicit in the previous theorem, but both 
metrics are conformal, with a conformal factor depending on a . By direct computation it follows 
that the curvature of this (Levi-Civita) metric is constant and equal to a. The checking for the 
last property could be even bypassed by suitably adapting and extending (so as to include also the 
parameter n 2 ) the idea used by Milnor for the Kepler problem in Euclidean space: consider the 
'inverted' momenta vector W, which may be called the 'slowmentum' and is defined as: 

Wi - 5f£b? m = sfrbt (50) 

Straightforward computations lead to several simple relations, whose standard k 2 = 1 versions are 
well known: 

W 2 + . 2 W 2 = 1 dW 2 + K2 dW 2 = ^ + K2 ^ 2 2 (51) 

allowing to obtain the expression of the Levi-Civita metric in terms of the 'slowmentum': 

4 iW )- t «,W) 
(l + <r(Wf + (aWJ)) 

which is precisely the form for the metric in a CK space with constant curvature a and signature 
type determined by k 2 in Riemann's normal coordinates. 

At this point, we expect that the well known connection between the Kepler motion in Euclidean 
space and the geodesic flow on constant curvature spaces extend so as to include also the Kepler 
problem on a 'curved' configuration space. This connection can be established following two paths 
which are however very closely related: either disclosing a stereographic projection implicit in the 
expressions we have found, or by a direct analysis of the geometry of cycles in the P\ , V-i-plane 
(compare Milnor [53] and Anosov [T] in the Euclidean case). In the next section we complete the 
details. 



3.1 Kepler motion as stereographic projection of free motion 

We start from the closed expressions for the 'momentum hodograph' corresponding to the CK orbit 
in the standard position. These follow by integration of the system (|44l) and with the correct choice 
of the integration constants they are: 

Vi{4>)=~ S K2 (0) , V 2 {4>) = (e + C K2 (0)) . (53) 
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This curve is a 'cycle' in the (Vi, T-^-plane, whose 'radius' is fc/( v / K2 l /) and whose center lies on 
the Vi axis, at the point (Vi, V2) — (0, fee/ (n-zJ)). This also determines the values of the constants 
foi,£'o2 for a Kepler orbit in standard position: 

c r, c k(l + e) , s 

£01 = 0, £02 = — • (54) 

The s-dependence of the CK momenta V\ , Vi can be found by replacing ([55]) in (|5"5|) : 

Vl{s) = -i-ec„( s y P2(s) = Ii- eCl ( s )' (55) 



where by direct computation: 



^ 2 W + ^2 2 W = ^ ? + e pf ! ■ (56) 
1 - e C CT (s) 



For the 'inverted' momentum or slowmentum vector, the s-evolution is: 

1 + e C CT (s) fc 1 + e C CT (s) 

and either by direct computation or using ([ST]) : 

Wl( S ) + K2 Wl(s)= l - \- eC /"\ (58) 
cr 1 + e C CT (s) 

The equations (|57[) can be directly recognized as the stereographic projection of a free motion in a 
space S^ K2 j with constant curvature a and metric of signature type K2 onto a space with a flat metric 
of signature type k 2 ■ Checking this statement involves a extension of the standard stereographic 
projection (where K2 = 1) to a /i2-general situation, making thus sense for all CK spaces. 

The classical stereographic projection maps the standard sphere S 2 of constant curvature 1 
(realized as the submanifold (s ) 2 + (s 1 ) 2 + (s 2 ) 2 = 1 in an auxiliar ambient space (s°, s 1 , s 2 )) on 
the flat plane s° = 1, s = w , s = w 2 living in this auxiliar space, by projecting the sphere from the 
'South pole' (s°, s 1 , s 2 ) = (—1, 0, 0); this mapping is well known to preserve angles. A geometrically 
similar construction maps any CK space S K i K2 \, realized as the 'sphere' (s°) 2 + ki(s 1 ) 2 + kiK2(s 2 ) 2 = 
1 in an auxiliar ambient space (s , s , s 2 ), over the flat plane s° = 1, s = w 1 , s 2 = w 2 , by projecting 
the 'sphere' over the plane from the 'South pole' (s°, s 1 , s 2 ) = (—1, 0, 0). This general stereographic 
projection is described, in terms of the ambient coordinates in S 2 by the map: 



( w 1 = 



l + s° 
s 2 

w = 



(59) 



1 + 8 







an expression which do not depends explicitly on k±,K2. In this form, the stereographic projection 
makes sense for all CK spaces, either with Riemannian or Lorentzian signature type. For H 2 , 
realized in the Weierstrass ambient space model, this stereographic projection provides the Poincare 
conformal disc model of the hyperbolic plane. 

Consider now in the CK space S 2 j K2 j with curvature a and metric of signature type K2, the 
'fiducial' geodesic described in ambient space coordinates as: 

(«W) = (C CT ( S ),-S ff ( S ),0) (60) 

which is the 'basic' line l\ along the direction 1 traversed negatively starting from the origin point 
(1, 0, 0). Then move this geodesic to a new position l[ by a translation along the orthogonal line I2 
with an amount e; we recall that e has a label UK2, just as the analogous quantity v has label K1K2 
in S K1 r K2 i . This translation map the fiducial geodesic l\ into some other member l[ of the set of 
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geodesies orthogonal to the basic line li- This one-dimensional family of geodesies is described in 
the ambient space by letting a general translation along the li line to act on the fiducial geodesic 
by matrix multiplication: 

C CTK2 (e) -™ 2 S CTK2 (e) \ / C CT ( S ) \ / C aK , (e) C^s) \ 

1 - Sa{s) = - Sa(s) (61) 

S CTK2 (e) C aK2 (e) J V / \ S™ 2 (e) C^s) / 
and this geodesic l[ in S^^] is mapped, upon stereographic projection (|59[) into the curve: 

l + C^te)^)' l + C CTK2 (e)C ff ( S ) l ° J 

to be compared with ((57j); these coincide provided the two identification conditions: 

Ca K2 (e)^e, S CTB2 (e) <-> ^ = )J - ^^k^ ^ 
are consistent. This requires 

CL 2 (e)+a K2 SL 2 (e) = l (64) 

and by direct checking we find this is satisfied, because this is simply (|13|) in disguise. The position 
of any geodesic in S^ri can be described by two parameters: an orientation angle <fio an d an 
'impact parameter' relative to the origin. The angle turns out to be precisely equal to the angle 
0o determining the orientation of the general Kepler conic in the configuration space. The 'impact 
parameter' e will contain, simultaneously, information on the eccentricity of the Kepler orbit and on 
the angular momentum (for the fixed value of the constant a). In particular, geodesies with impact 
parameter equal to zero correspond to collision orbits with J = 0; from the geometric picture, as 
the symmetry group of SjL^n acts transitively on points and on either type ('time' or 'space' like 
when K2 < 0) geodesies, it follows that the collision orbits receive, alike any other orbit, a regular 
description within this scheme. 

The picture is the following: On a CK configuration space S K i K2 \ with constant curvature K\ 
and metric of signature type K2, the subset of all Kepler motions for which the quantity a = 
— (2E — K1K2J 2 ) has a fixed value can be realized as the geodesic flow on a CK space S^ K2 j of 
constant curvature a and metric of the signature type Ki. When a (fictitious) point on this space 
moves with unit speed along a geodesic l[ obtained from the 'basic' geodesic l\ by a translation 
of amount e along I2, its stereographic projection on the 'inverted' momentum space corresponds 
to a (slowmentum hodograph) of the Kepler motion along a orbit in the standard position, with 
eccentricity, semilatus rectum and angular momentum determined from e by (|63p . The arc lenght 
along the geodesic l[ is the 'curved' Levi-Civita parameter. 

Let us rewrite the results of the previous section, in terms of the identification (f6"3"| : 

(T Kl H)(.s) = fc Cj(s) ~ CgK2(£) , (T KlK2 ( W ))( S ) = fcS CTK2 (e)S (T ( S ), (65) 

(T Bl (r))(«) - fc i-eWOfirfr) (tm){s) = fcS CTK2 (e) S g ( a ) a r (66) 

The structure of these expressions is neat, and their genericity must be emphasized: they hold 
for any Kepler motion, with any a, in any S 2 ^ K2 ^ . All involve as a factor the strenght of the 
Kepler coupling constant fc. In all cases there are two further independent variables, which can be 
identified to parallel T' coordinates in SV,. One of these coordinates corresponds to the choice 
of a standard geodesic (the 'impact parameter' e) and the other is the arc lenght parameter along 
the geodesic. The presence of explicit a in the denominators would seem to imply some difficulties 
when a —> 0, but this is not so. As a consequence of the behaviour of the Cosine and Sine functions, 
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the vanishing of a is automatically accompanied by the vanishing of 1 — C a (s). By introducing the 
versed sine V CT (s), defined for nonvanishing a as: 

Vff ( S ) = l^£l, (67) 

cr 

then when cr — ► this function reduce to s 2 /2 and should be considered as the next natural 
stages in the 'curved' analogues of s 2 /2 just as C CT (s) and Ser(s) are the 'curved' analogues of the 
functions 1 and s. When (|651 I66p are reexpressed in terms of the functions V CT (s) and V CTK2 (e) we 
obtain a description of the s-evolution which covers the Kepler motion on any CK space (with any 
constant curvature and any signature type) and which makes furthermore sense for any value of 
the combination a of energy and angular momentum: 

(t ki (u))( S ) -fc( K2 V CTK2 (e)- V CT (s)) , (T Kire »)(s) = fcS CTK2 (e)S ff (s), (68) 
(T Kl (r))( S ) =fc( K2 V. K2 ( £ ) + C CTK2 (e)V ff ( S )) , ( T^)) (s) = fc ^'{ff (69) 

v / v y v y K2V CTK2 (e) - v CT (s) 

In particular, orbits with a = 0, which are in some aspects analogues in the 'curved' configuration 
space of the (non-generic) parabolic orbits in Euclidean space, are described as the particular case 
cr = of the previous generic expressions, which simplify in this case to: 

( T Kl (u)) (s) = k - j V ( T KlK2 (v)) 00 = fees , ( T Kl (r)) (s) = k + ^ , 

(70) 

so, for these orbits with a = and no matter of the value of the curvature of the configuration 
space, the s-evolution is the following: T KlK2 (v) is a linear function of s and T Kl (u) and T Kl (r) are 
quadratic (further t is cubic). This extends the semicubical type singularity at the cusp for the 
function x{t), whose graph is a cycloid for the parabolic orbits in the Euclidean case. 

A point worth to remark is the following: for any a and in a configuration space of any curvature, 
T KlK2 (w) is still a linear function of S CT (s), which is itself the 'cr-deformation' of s, so in some CK 
sense l T KlK2 (v) is linear in s', while T Kl (u) and T Kl (r) are linear functions of V CT (s) which the is the 
natural 'cr-deformation' of s 2 /2, hence they are 'quadratic in s' in the CK sense. 

3.2 Kepler motion on a curved space 5 , ^j K2 ] as a geodesic flow in a S^ K2 j 

The results obtained in the previous sections can be summed up in the form: 

Theorem ['curved' version of Moser, Osipov, Belbruno; Milnor] Consider Kepler motions of a 
particle in a configuration space S^^ K2 ^ with a metric of constant curvature K\ and signature type 
K2 (Riemannian or Lorentzian for k 2 >, < respectively). For such a Kepler motion, there are 
two basic constants of motion, the energy E and the angular momentum J . Restrict attention 
to the set of all Kepler motions with a fixed value of the combination a = —(2E — K1K2J 2 ), and 
look to these motions in momentum plane V — {Vi,V2)- On the space of all 'momentum vectors' 
admissible for such motions, suitably completed , there is a unique metric ds 2 , which is of signature 
type K2 , with the following properties: 

1 ) The metric has constant curvature a. 

2) The geodesies of this metric are precisely the 'momentum hodographs' for the 'curved' Kepler 
problem, and 

3) The arc-lenght parameter of this metric coincides with the 'curved' Levi-Civita parameter 
along the corresponding Kepler motion. 

Hence, the Kepler motion on a curved space S 2 _^ K2 ^ can be seen as a geodesic flow in a S 2 r j 

Two comments would help to clear up possible misunderstandings. First, the 'suitable comple- 
tion' mentioned in the Theorem statement is the following: When k 2 > 0, a single point at infinity 
V = 00 has to be added. When k 2 < 0, the completion requires to add a single point at infinity 
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and two straight lines, which together make the light cone of the point at infinity; some detailed 
discussion of this issue in [6;. Second, the previous theorem refers motions with a fixed given a 
and to a metric, the 'Levi-Civita one' ds 2 which is of constant curvature cr, in the (completed) 
momentum plane V = (Vi,V2), turning this space into a S^ Ka j. In the same momentum plane, 
without the completion, we may consider another flat metric with signature type K2 as well and 
independent of cr. The momentum hodographs of any Kepler orbit, which are geodesies relatively 
to the Levi-Civita metric with the corresponding value for cr, can be alternatively seen as 'cycles' 
relatively to the flat metric dp 2 = V\ + K1V 2 . For a fixed value of cr, the cycles are precisely the 
ones obtained from the stereographic projection of the geodesies in S^, re 1. Of course, completion 
of this plane is required if we want to have an uniform description including e.g., straight lines as 
circles. 



3.3 The period of the 'curved' Kepler orbits and the three Kepler laws 

The major semiaxis a — i(r per + r apo ) of the closed elliptic orbits (those with e < 1 and real r per 
and r apo ; K2 > 0) can be easily related to the energy and angular momentum of the orbit. We start 
from (jTSJ) , evaluate it at 4> = (the periastron) and at cj> = ^/\f^2 (the apoastron), and use (fTTIf 
to obtain: 

T ( \ — K2 ^ 1 T I \ — K2 ^~ 1 (T\\ 

J- K^per ) — ^ 1 + e' L Ki\ r apo) — ^ 1 — e' 

Then expand T Kl (2a) = T Kl (r per + r apo ) using the formula for addition of tangents (|15j) 

rn / , \ T Kl (ri) + T Kl (r 2 ) 

1 - Ki T Kl (ri) T Kl (r 2 ) 

and the relation (|13p . Angular momentum dissapears and the result is 

T K1 (2a) = (73) 

hence for the closed elliptic orbits which appear in Riemannian configuration spaces of any curvature 
«i, the energy of the motion depends only on the ellipse major semiaxis, and the relation is: 

T Kl (2a) 

which evidently reduces in the Euclidean case k% = to the well known result E — — k/(2a). 

The period for these orbits can be also similarly obtained. Start from noticing that periastron 
and apoastron correspond to the two values s — and s = 2 \^_ a = n/^/a, so for symmetry reasons 
the period of a closed elliptic orbit can be directly obtained as T = 2t{s = 2\^ a ) in terms of formula 

(|4T|) . Now, as T CT (^ CT ) = oo, multiplication by j£\ is ineffective, and ArcT_ CT ^-^=±-L= T CT (^ c 

[i. . The same happens for the (3 contribution. Thus we get the relation 

ekni T= J a \ ir 

2a \Vc^T ^~^\\2^' { ' 

Now if we develop this expression using (|39p and (|7^)) . after some slightly tedious but straight- 
fowward computation, the angular momentum also disappears and the period depends only on the 
energy, as was to be expected, through a rather uninspiring relation: 



T 2 = n 2 n- ^ (76) 



and in the standard spherical case this coincides, after some algebra, with the one given in [18] . 
Again when k\ — > the limit is easily seen to be T 2 = tt 2 z!km ■ Back to the general case with 
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any curvature «i, by using the relation (|74p between energy and major semiaxis and simplifying, a 
much more transparent relation relating the major ellipse semiaxis to the period appears: 



7T 2 2lT 2 47T 2 

T 2 = — S Kl (2a) V Kl (2a) = — S Kl (2a) S 2 Kl (a) = — C Kl (a) S 3 Kl (a) , (77) 

which reduces in the flat Euclidean configuration space to the classical result 

2 , ir 2 (2a) 2 4Tr 2 a 3 , s 

For the standard hyperbolic space, the last expression in ([77]) was obtained also by a direct com- 
putation starting from the law of areas by Liebmann [22] (in the chapter entitled "Nichteuklidische 
Mechanik" of the 1905 edition). Probably it is more clear to stick to the first form, as it displays 
a factor 'CK linear' and a factor 'CK quadratic' in the same variable 2a, in the sense discussed 
before. In terms of the natural pulsation lu := 2w/T for the Kepler orbits, these relations can be 
written, for any CK space, resembling the Euclidean form of the third Kepler law known as the 
1-2-3 relation: 

k = C Kl (a)u 2 S 3 Kl (a) (79) 

(the term C Kl (a) becomes invisible in the flat k% = case, where the formula reduces to k = u> 2 a 3 ). 

In all the previous sections, most Euclidean expressions involving E have a 'curved' analogous 
which is formally similar when expressed in terms of a instead. The relations between energy and 
ellipse major semiaxis or period do not follow this pattern , as their 'curved' forms involve directly 
the energy E, just as in the flat case. 

Hence we may formulate the three classical laws for the Kepler problem on any configuration 
space S Kl r K2 i in the following form: 

• Kepler orbits are conies, with a focus at the potential center. 

• The law of areas, in the form J = 8 Kl (r)(j> is constant, always holds. 

• For closed elliptic orbits, period T (or pulsation u>) and semimajor axis a are related by ([771 179")) . 



4 From the 'curved' to the Euclidean Kepler problem 

4.1 The Kepler problem on the Euclidean space 

Of course, setting the standard values n\ = 0, K2 = 1 for the Euclidean space, all the expressions in 
Section 2 reduce to the corresponding ones for the Euclidean Kepler problem. All these expressions 
are quite well known, and the only novelty remaining in the Euclidean case lies in the use of a 
unified setting, valid for any value of the energy, by means of the parameter a instead of the energy 
E = —a/2. Indeed most of the general 'curved' expressions are very similar, mutatis mutandis, 
to the Euclidean ones, and the only difference is the general appearance of labelled trigonometric 
functions of either u, r, y, <j> with its natural labels, respectively k\, k\, k\K2, hi which reduce in the 
Euclidean case to parabolic functions of u = x, r, y and to circular functions of </>. For instance, for 
the relation giving the s evolution of u: 

T Kl (u) = ^ ( C a (s) - e) reduces in E 2 to u(= x) = ^ ( C^s) - e) . (80) 

The s-dependence of the cartesian coordinate y, the radial coordinate r and of the angular coordi- 
nate <f> similarly reduce to: 

y(s) = JS a (s), r{a) = ~ex(a) = ^(l-eC^a)), tan = (1+e)^-^ T a (|) . (81) 

The only not directly obvious Euclidean limiting form is the s-dependence of the time t(s), 
where a cursory look to (|4"Tj) does not disclose resemblance to any familiar Euclidean formula. The 
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naive replacement k\ = in (141|) gives an indeterminate expression of the type 0/0. The reason 
of course comes from the K\ dependence of the integrand in ()37p : when n\ — the integration is 
immediate and does not call for a decomposition into simple fractions. A more careful analysis of 
the situation involves dealing with the K\ — > limit. This is an exercise which requires a careful 
consideration of square root determinations and signs, but is otherwise straightforward, so details 
will not be given here. The result is the one one should expect: in terms of a power expansion in 
K i, (EH) appears as: 

= k s-eS„(s) + 
a 

which directly recalls the Kepler equation. The precise connection will be discussed in the next 
section. 



4.2 The Levi-Civita parameter s versus the eccentric anomaly 

Traditionally, the parameter used to describe the evolution along Euclidean Kepler orbits, the 
eccentric anomaly, denoted £, can be introduced for either type of orbits (here we follow [28| 12]) 
as: 

x = a(cos£ — e), y = a a/1 — e 2 sin£ when E < 0; a = -y^ 

x=i(p-?), V=y/§Z, when^-0 (83) 

x = a(cosh£; — e), y = as/l — e 2 sinh£ when E > 0; a = ^ . 

The parameter £ is defined so that, for any negative energy, a complete revolution along the orbit 
corresponds to an increase of £ from to 2ir. When the energy is positive, £ takes all real values, 
and its scaling is chosen in analogy with the case of elliptic orbits. By replacing x(s), y(s) as taken 
from (|80l l8Tj) and comparing, one finds that s and £ are related by simple scalings: 

£ = a/cxs, cos£ = C CT (s), sin£ = y/trSJ^s) when E < , 

£ = v/fes, when E = , (84) 

£ = a/— its, cosh£ = C(j(s), sinh£ = a/— <t S^s) when i? > , 



so that along any orbit, the eccentric anomaly £ is proportional to s, but the coefficient depends on 
the value of the orbit energy, and is not continuous at E = (even dimensionally this 'parabolic' 
eccentric anomaly differs from the 'elliptic' and 'hyperbolic' ones and furthermore, different authors 
define it differently, compare [19j). All the 'universal' equations in the previous section involving s 
can be rewritten in terms of £, but this conceals (specially when E — 0) the underlying similarity. 
From our present viewpoint, this can be seen as an unwanted artifact of the normalization imposed 
to the eccentric anomaly for the general elliptic and hyperbolic orbits, in a way which precludes a 
non-trivial limit for £ when a — > (as clear in l|84p). For instance, the Euclidean equations for t(s) 
([82l with Ki = 0) and </>(s) (|8T|) reduce in three three types of orbits with E <, =, > respectively 
to: 




tan- 



'-(esinh£-£), 



1 



1 



— tan — , 

e 2 



3 



(85) 



^-Uanh|, 
e — 1 2 



Summing up: the facility of using the unsealed circular or trigonometric functions of the eccentric 
anomalies in the description of all orbits with non-zero energy, as made evident in ()83l I85p , comes at 
the price of losing the manifest universality of the approach using s as a basic variable instead. In 
this approach, the functions of s which appear naturally have to be the scaled associated functions 
C CT (s), S a (s) whose CK label is a = —2E but this description remains smooth when E — > 0. 
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